Absence of Wavepacket Diffusion in Disordered Nonlinear Systems 
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We study the spreading of an initially localized wavepacket in two nonlinear chains (discrete 
nonlinear Schrodinger and quartic Klein-Gordon) with disorder. Previous studies suggest that there 
are many initial conditions such that the second moment of the norm and energy density distributions 
diverge as a function of time. We find that the participation number of a wavepacket does not diverge 
simultaneously. We prove this result analytically for norm-conserving models and strong enough 
nonlincarity. After long times the dynamical state consists of a distribution of nondecaying yet 
interacting normal modes. The Fourier spectrum shows quasiperiodic dynamics. Assuming this 
result holds for any initially localized wavepacket, a limit profile for the norm/energy distribution 
with infinite second moment should exist in all cases which rules out the possibility of slow energy 
diffusion (subdiffusion). This limit profile could be a quasiperiodic solution (KAM torus). 
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PACS numbers: Valid PACS appear here 



It is well-known that Anderson localization occurs for a 
one-dimensional linear system with uncorrelated random 
potential. Since all the linear eigenmodes - Anderson 
modes (AMs) - are localized, any wavepacket which is 
initially localized remains localized for all time. There- 
fore there is no energy diffusion When nonlineari- 
ties are added to such models, AMs interact with each 
other, giving rise to more complex situations Q • Numer- 
ical studies of wavepacket propagation in several models 
showed that the second moment of the norm/energy dis- 
tribution growths subdiffusively in time as t a [1, 0, @, 
with a in the range 0.3 — 0.4, though not being accu- 
rately determined. The conclusion was that the initial 
excitation will completely delocalizc for infinite times. 
Recently, experiments were performed on light propaga- 
tion in spatially random nonlinear optical media [6| . 

Spatially periodic nonlinear systems will support Dis- 
crete Breathers (DBs), which are spatially localized time 
periodic solutions Q with frequencies outside the fre- 
quency spectrum of the linear system. The temporal 
evolution of a localized wavepacket leads to the forma- 
tion of a DB, while a part of the energy of the wavepacket 
is radiated ballistically to infinity (in the form of weakly 
nonlinear plane waves) Q. In that case, the second mo- 
ment of the energy density distribution diverges as t 2 , 
falsly suggesting complete derealization. The participa- 
tion number P of the norm/energy distribution (or sim- 
ilar quantities) is a well-known measure of the degree of 
localization. In the case of a periodic nonlinear lattice, 
P will saturate at a finite value, correctly indicating the 
formation of a DB. 

For nonlinear random systems it was proven rigorously 
that AMs survive in the presence of nonlinearities as spa- 
tially localized and time-periodic solutions @ with fre- 
quencies which depend on the amplitude of the mode. 
The allowed frequencies form a fat Cantor set (with finite 
measure) whose density becomes unity for weak nonlin- 



carity. They are located inside the frequency spectrum 
of the linear system. Numerical techniques for obtain- 
ing these (dynamically stable) intraband DB solutions at 
computer accuracy were developed [lOt ] . When they are 
chosen as an initial wavepacket, they persist for infinite 
time and there is no diffusion at all. 

Here we analyse carefully the evolution of the partic- 
ipation number of wavepackcts as a function of time, in 
situations where subdiffusion is claimed to exist [!, 0, [f| . 
We study two models. The Hamiltonian of the disordered 
discrete nonlinear Schrodinger equation (DNLS) 
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with complex variables ip n . The random on-site energies 
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e n are chosen uniformly from the interval In- 
equations of motion are generated by ip n = dHD/d(iip*). 
We choose j3 = 1 and V = — 1 here [J| and note that vary- 
ing the norm of the initial wavepacket is strictly equiva- 
lent to varying (3. 

The Hamiltonian of the quartic Klein-Gordon chain 
(KG) 



H * = E Y + f < + \9< + j(Un+l - Unf . (2) 
n 

The equations of motion are «„ = —dH.K /du ni e n — 
1 + e n (W = 1), and g = 1. 

For (3 = g = both models are reduced to the linear 
eigenvalue problem XA n = e n A n — V(A n+ i +A n -\). The 
eigenvectors A v n are the AMs, and the eigenvalues A„ are 
the frequencies of the AMs for the DNLS, while the KG 
modes have frequencies oj v = VA„ + 1 + 2V. 

Hamiltonian |T|) (unlike @), in addition to conserving 
the energy, also conserves the total norm S = J2 n IV'nl 2 = 
(ip\ip). We use this norm conservation for proving rig- 
orously that initially localized wavepackets with a large 
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enough amplitude cannot spread to arbitrarily small am- 
plitudes. The consequence is that a part of the initial 
energy must remain well-focused at all times. 

This proof is inspired by [ll[ . We split the total energy 
Hd = (V'l-MV') +-^iVL into the sum of its quadratic term 
of order 2 and its nonlinear terms of order strictly higher 
than 2. Then, L is a linear operator which is bounded 
from below (and above). In our specific example, we 
have (t/>|L|i/>) > £ * ; m(V'lV') = ^mS where cu m > —2 — ^ is 
the lowest eigenvalue of L. Otherwise, the higher order 
nonlinear terms have to be strictly negative. 

If we assume that the wavepacket ampli- 
tudes spread to zero at infinite time, we have 
lim t _ >+oo (sup n |'0n|) = 0. Then lim t _ >+o0 (X)„ |"0n| 4 ) < 
lim^ + oo(sup„|^i)(E„IV'n| 2 ) = since S = £ n V£ 
is time invariant. Consequently , for t — > +oo we have 
Hnl = and Hd > u rn J2 n \^"\ 2 = u m S. Since Hd 
and S are both time invariant, this inequality should be 
fulfilled at all times. However when the initial amplitude 
A of the wavepacket is large enough, it cannot be 
initially fulfilled because the nonlinear energy diverges 
as —A A while the total norm diverges as A 2 only. For 
example, a wavepacket initially at a site ( tp n = for 
n / and ipQ = \f~A ) has energy Hd = £qA 2 — \A A . 
Consequently, the above inequality is not fulfilled when 
A 2 > —2(u m — cq) > 0. Thus such an initial wavepacket 
cannot spread to zero amplitudes at infinite time. 

This proof is valid for DNLS models with any W (in- 
cluding the periodic case) and any lattice dimension and 
can be easily extended to larger classes of DNLS models 
where the nonlinear terms are either strictly negative, or 
strictly positive. Note that the large amplitude regime 
where we prove that complete energy diffusion is impos- 
sible in DNLS models, is precisely the one where subdif- 
fusion is claimed to completely delocalize the wavepacket 
P|. Thus we disprove these claims. 

We performed extensive numerical simulations, and 
characterized the wavepacket spreading both in real 
space for DNLS and normal mode space (Anderson space 
or AS) for KG. We used initial wavepackets with all the 
energy localized on a single site no, or single AM, or com- 
binations, close to no. Nonlinearity induces diffusion in 
Anderson space, where each AM is characterized by a am- 
plitude a u and momentum a„. We analyze distributions 
zi > using the second moment mi = — Iq) 2 zi and 
the participation number P = (J^ ; zi) 2 j z 2 , which 
measures the number of the strongest excited sites in 
Zi. We order the AMs in space by increasing value of the 
center-of-norm coordinate X v = J2 n nj ^n- I n the results 
presented here, for the DNLS z n — \ip n \ 2 is the norm den- 
sity in real space, and for the KG z v — a 2 /2 + tu 2 a 2 /2 
is the (harmonic) energy density in AS. The system size 
was N = 1000 for KG, and N = 2000 for DNLS. Excita- 
tions did not reach the boundaries during the integration 
time, and results are unchanged when further increasing 
N. 

We show in FigfT] the KG energy distribution in AS 
for a single site excitation with energy E — 1, and 
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FIG. 1: (color online) KG: Energy distribution at t = 6 x 10 7 
(black solid) and t = 1.2 x 10 8 (red dashed) in AS. Initial 
single site excitation with energy E = 1, V = 0.25. Insets: 
profiles of the strongest excited AMs in real space. 



V = 0.25 at times t = 6 x 10 7 , 1.2 x 10 8 . Two rather 
strongly excited modes are surviving almost unchanged 
on these time scales. The insets show their eigenvectors, 
which are well localized, and practically do not overlap. 
The same distributions on a logarithmic scale (KG and 
DNLS) show a chapeau of weaker excited AMs, with ex- 
ponential tails due to its finite width (Figf2|). This cha- 
peau is perhaps slowly growing. The subdiffusive growth 
of the second moment at these times (see FigJ3]) is mainly 
due to weak excitation of tail modes. 

The participation number P(t) is plotted in Figl3]for 
the same runs. We observe no growth. P fluctuates 
around a value of 7-10, confirming the results in FigfTJ 
that we observe a localized state, similar to a DB. As- 
sume that the rest of the weakly excited modes continues 
to subdiffuse in the chapeau. We use a modified distribu- 
tion z v for the KG run, where the 10 strongest mode con- 
tributions are zeroed (top panel, green curve). The weak 
mode participation number is now fluctuating around 70, 
but again does not grow. Therefore the chapeau appears 
not to diffuse, and the observed growth of 7712 ~ £0.3.. .0.4 
is not related to a derealization process. Instead, we 
find that the packet does not delocalize. Indeed, assum- 
ing that the chapeau homogeneously spreads in a sub- 
diffusive way as claimed, it follows that P(t) ~ i Q / 2 , 
which clearly contradicts our observations. We repeated 
these runs with various initial conditions and disorder re- 
alizations with similar results. However the localization 
pattern (Figfl]), and the observed averaged participation 
number P, fluctuate. Performing an averaging of the final 
distribution over several realizations 0, U IM w * n there- 
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FIG. 3: (color online) P and m.2 versus time, on logarithmic 
scale. Parameters as in Fig[2] Top panel: KG, AS. Bottom 
panel: DNLS, real space. 

fore completely smear out the sharp localization patterns 
in the distributions. Closer inspection of the evolution of 
m2 shows, that the exponent a is strongly depending on 
the time intervals of study, and also on the given disorder 
realization. There are some indications suggesting that a 
might decay at long time and even that ui2 (t) may satu- 
rate, but further clarification may call for very extensive 



numerical investigations. 

Finally we calculated the Fourier transform I(uj) of 
P(t) (after ( = 2x 10 7 , over an interval of At = 2000), 
see FiglU We find a quasiperiodic spectrum, which is 
close to periodic, with no hints of a chaos-induced con- 
tinuous part. For the KG case the energy densities are 
quadratic forms of the AM coordinates, thus the main 
peak position ui « 3 corresponds to a frequency u) « 1.5 
for the AM coordinate dependence, which coincides with 
the frequencies of the strong excited Anderson modes in 

Figru 

Our main result is, that in both models |T]) and ([2]), 
whatever the initial wavepacket is (even if it is not fulfill- 
ing the conditions for our theorem), and irrespective of 
the model parameters and the disorder realization, the 
participation number does not diverge as a function of 
time as it should in case of subdiffusion (as t a / 2 ) but in- 
stead fluctuates between finite upper and lower bounds. 

Let us now propose an interpretation of our observa- 
tions. First, it is useful to recall the wavepacket behavior 
in the absence of disorder. When its amplitude is large 
enough for generating a DB, there is a transient dynam- 
ical state which is more or less chaotic, with a broad 
band time-Fourier spectrum overlapping the spectrum of 
the linear system. Because of that, a part of the energy 
of the wavepacket is radiated to infinity. With that, the 
remaining DB like excitation becomes quasiperiodic first, 
and finally, approaches an equidistant spectrum of peri- 
odic motion, which completely stops further radiation. 
The energy which has been emitted spreads towards in- 
finity. Therefore there is a limit profile which is a lo- 
calized time-periodic solution - an exact DB. This is the 
only possibility for the limit profile, in order to avoid ra- 
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diation. This is an example where the initial wavepacket 
self-organizes in order to stop radiation. 

When the system is both random and nonlinear, ra- 
diation into the linear spectrum is impossible due to 
Anderson localization. Nevertheless, the same process 
starts as before, but the energy emitted by the initial 
wavepacket cannot spread towards infinity since the par- 
ticipation number (full and partial) does not diverge. 
The following cascading scenario may be true. The core 
of the wavepacket emits a part of its energy which re- 
mains within the linear localization length nearby the 
initial wavepacket (due to the nonlinearity-induced cou- 
pling between the AMs). The same process should re- 
peat for the emitted energy. A part of it remains local- 
ized while another part is reemitted a bit farther from 
the central site within the localization length and so on. 
This process of reemission repeats forever and generates 
a tail for the wavepacket which will become much more 
extended than the localization length. The central am- 
plitude of the wavepacket does not tend to zero. The 
process of energy reemission slows down when the ampli- 
tude at the edge of the tail becomes small which explains 
the very slow numerical convergence. The final result is 
that at infinite time, the energy (or norm) distribution 
should converge to a nonvanishing limit profile which is 
summable since energy (or norm) is conserved. However, 
it may or may not have a finite second moment, which 
makes the question of the evolution of the second moment 
secondary. Unlike the standard DB case in spatially pe- 
riodic systems, the limit profile is not a time periodic 
solution. 

It was proven rigorously ([H,[l3j]) that stable spatially 
localized quasiperiodic solutions with finite energy ex- 
ist in similar nonlinear models with infinitely many de- 
grees of freedom without or with degenerate linear spec- 
trum. These KAM tori are quasiperiodic DBs which in 
some sense are linear combination of Anderson modes 
surviving in the presence of nonlinearity. Indeed, we find 
that the Fourier spectrum of the wavepacket dynamics 
becomes quasiperiodic, with narrow peaks and a small 
background as time grows suggesting the motion tends 
to become quasiperiodic (FigHJ). 

If the limit profile becomes a KAM torus, we should 
also observe that the largest Lyapunov exponent tends 
to as t — » +00. Indeed, we find that this Lyapunov 



exponent drops rapidly during the first expansion part of 
the wavepacket, and slowly further decays, with charac- 
teristic values of 10~ 4 at the end of our simulations. The 
corresponding time scale is 10 4 , and four orders of mag- 
nitude smaller than the simulation times. No chaotic 
dynamics is observable, and we think that the conver- 
gence to the final KAM torus is very slow because the 
surrounding KAM tori are expected to become dense. 
We should even expect to enter the regime of Arnol'd 
diffusion which is expected to be very slow and difficult 
to investigate both numerically and analytically. 

Note that this convergence to a quasiperiodic limit pro- 
file can only occur in infinite systems because if the sys- 
tem is finite, the regularization process of the initially 
chaotic trajectories ends when the packet tails reach the 
edge of the box. Then, we should expect to get equipar- 
tition of the energy after a sufficiently large time and a 
trajectory which remains chaotic with a nonzero largest 
Lyapunov exponent. 

In summary, we have proved by a rigorous analytical 
argument, and completed by numerical investigations of 
the participation number, that a wavepacket in a ran- 
dom nonlinear system does not spread ad infinitum. A 
limiting quasiperiodic profile is approached, and the slow 
increase of the second moment of the energy/norm dis- 
tribution does not violate these findings. It is an open 
question whether the limiting profile will have a finite 
or infinite second moment. Thus, we observe absence of 
diffusion in nonlinear disordered systems. Note that this 
conclusion can be equally well applied to higher dimen- 
sional systems, provided all AMs are localized. 
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